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Abstract 

The gauged Nambu- Jona-Lasinio model in the quenched-ladder approxima- 
tion has nontrivial dynamics near a critical scaling region (critical curve) 
separating a chiral symmetric and a dynamically chiral symmetry broken 
phase. Scalar and pseudoscalar composites corresponding to the four-fermion 
and gauge interactions become relevant degrees of freedom at long distances, 
which is reflected in the appearance of a large anomalous dimension of the 
four-fermion operators. A method is introduced for solving the Schwinger- 
Dyson equation for the Yukawa vertex in specific kinematic regimes. This 
allows one to derive an analytic expression for the scalar propagator, which is 
valid along the entire critical curve. The mass and width of the scalar com- 
posites in the critical scaling region are reexamined and the conformal phase 
transition at gauge coupling qq = Oc is discussed. 

ll.lO.St, 11.15.Tk, ll.30.Qc, 11.30.Rd 
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Typeset using REVT^ 



I. INTRODUCTION 



The recent interest in the gauged Nambu-Jona-Lasinio (GNJL) model has been stimu- 
lated by its importance for constructing extended technicolor models (ETC) and the top- 
quark condensate model (for an introduction see Ref. 

It is well known that in the quenched-ladder approximation]^ the GNJL model has a non- 
trivial phase structure in the coupling constant plane {go, ao), where go = GqA'^/Aii^, = 
el/An as shown in Fig. |1| (A is a cutoff). The critical line is 




< «o < ttc = ^ (1) 



at go > 2, and 



ao = ac (2) 

at (7o < |) above which the gap equation for the fermion self-energy S(p) has a nontrivial 
solution. Thus the chiral symmetry is dynamically broken, which implies the existence of 
a nonzero vacuum condensate {ipip)- One end point {go = l,ao = 0) of the critical line 
corresponds to the ordinary NJL model, while the other one {go = 0, ao = «c = 7r/3) 
corresponds to pure QED. The interesting feature of the GNJL model is the observation 
that naively irrelevant chiral invariant four-fermion operators become relevant near the chiral 
phase transition due to the appearance of a large anomalous dimension 7^ of the operator 
%l)ip along the critical curve Eq. ([l|), 2 > 7^ > 1, and 7^ = 1 along the part of the critical 
line with ao = @]- Therefore the GNJL model is believed to be renormalizable as an 
interacting continuum theory near a critical scaling region (critical curve) in the coupling 
plane separating a chiral symmetric phase {xS) and a spontaneous chiral symmetry broken 
phase {SxSB). 

In the ladder approximation it has been shown that the GNJL model in four 

dimensions is indeed renormalizable, and that the anomalous dimension 

_ (91nmo , . 

7m = -A^^ (3) 

is 

7m = 1 - + 2a;— , 5-0 < fi'c = = -y/l - ao/ac, (4) 

7m = l + a; + 2(5(0 go>9c- (5) 

Fine-tuning the coupling g^ to gc in SxSB phase in such a way that rrid/A <^ 1, where 
rrid = S(0) is the dynamical mass of a fermion, a nontrivial continuum limit (m^/A — > 0) 
can be reached just as in pure quenched QED p-pHj]. The spectrum of such a theory contains 



pseudoscalar (vr) and scalar (cr) bound states which become light and dynamically active in 



^Also referred to as the quenched-planar approximation. 
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the vicinity of the critical hne. Since the phase transition is second order along the part 
Eq. dll) of the critical curve, scalar and pseudoscalar resonances have been shown to be 
produced on the symmetric side of the curve, whose masses approach zero as the critical 
curve is approached |T^. The part of the critical curve, Eq. (|]), with = etc is rather 
special. For example, an abrupt change of the spectrum of light excitations occurs when the 
line ao = Oc, (70 < 1/4 is crossed: while light scalar and pseudoscalar excitations still persist 
in the broken phase, there are no such light excitations in the symmetric phase P,p!^. A 
similar behavior has been revealed also in QED3 |T^. This peculiar phase transition was 
referred to conformal phase transition (CPT) | |14| . 

In this paper we study scalar composites {a and tt bosons) in the symmetric phase of the 
GNJL model. Computing the scalar propagator (see Fig. |]) requires knowledge of the full 
scalar-fermion-antifermion vertex ^s{p + (l,p) which in turn satisfies the Bethe-Salpeter (BS) 
equation displayed in Fig. ||. We solve this BS equation in the ladder approximation (Fig. H) 
but differ from the corresponding studies in Refs. 0JTT| who used an approximation for 
rs(p + g,p) with zero boson momentum (g = 0). A technique of expansion in the Chebyshev 
polynomials is introduced for solving the Yukawa vertex with nonzero boson momentum and 
consequently an explicit analytical expression is derived for the propagator of the a boson 
valid along the entire critical curve. 

Our main physical conclusions are the same as in Refs. |^,11,12|: in the region ao < Oc, 



go > 1/4, in the symmetric phase, a spectrum of light resonances exists while at ao < 
aciQo < 1/4 there are no light resonances. Having obtained an analytical expression for the 
scalar propagator, we can analytically continue it into the region > Oc and find light 
tachyons there, signaling the instability of the symmetric solution. 

The plan of the present paper is as follows. After introducing the GNJL model (Sec. H) we 
solve the equation for the Yukawa vertex with nonzero boson momentum in Sec. ^ keeping 
only the zero order Chebyshev harmonics. In Sec. |^ we obtain an analytical expression for 
the boson propagator valid along the entire critical line and analyze its behavior in different 
asymptotical regimes. Sec. is devoted to comparing our results for the Yukawa vertex 
and boson propagator with the corresponding ones in Refs. [l7|,|Tl|. In Sec. we discuss 



the behavior of the scalar propagator near the critical line ([^) in the symmetric phase and, 
in particular, the mass and the width of resonances. The analysis of the scalar composites 
near the critical line (^ is given in Sec. |VII| , where we show the absence of light excitations 
at tto < otc while analytically continuing the symmetric phase propagator into the region 
Q^o ^ otc leads to the appearance of tachyonic states. We discuss this behavior from the 



viewpoint of the CPT conception proposed recently in Ref. |TJ]. We present our summary 
in Sec. |V111 and in Appendix ^ we give an analysis of the contribution of higher order 
Chebyshev harmonics into the Yukawa vertex equation and scalar vacuum polarization. 



II. THE GAUGED NAMBU JONA-LASINIO MODEL 

The gauged NJL model is described by the Lagrangian 

C-GmL = -\fI, + ^{il^D^ - m,)^ + + {i^i^.^fl (6) 

where = — ie^A^ is a covariant derivative, Cq the gauge coupling constant, and the 
last term is a chirally invariant four-fermion interaction with Gq the corresponding Fermi 
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coupling constant. Another way to refer to this model is that it is QED with an additional 
four-fermion interaction. In the absence of a fermion mass term mg which breaks the chiral 
symmetry explicitly, the Lagrangian (||) possesses a U{1) gauge symmetry and a global 
^7^(1) X Uji{l) chiral symmetry. (It is not difficult to extend our results for SUii^Nf) x 
SUii{Nf) chiral symmetry with Nf fermion flavors.) 

Let us introduce the chiral fields a and vr rewriting the Lagrangian (^ in the form 

^GNJL = -]f^ + i^i^^D^i, - ij{a + ^757r)V^ - :^(^' + + 77^^- (7) 

One can readily verify the equivalence of the Lagrangians (P) and (|^) by just making use 
of the Euler-Lagrange equations. Further we study mainly the chiral symmetric case with 
niQ = 0. 

With this Lagrangian we define the generating functional of the GNJL model by 

2'gnjl[^, V, V, JcT, Jn] = N J d/i (^A, ijj, ijj, cr, vr) exp {iSj} , (8) 
where the "measure" is defined as 

d/i (a, ip, ip, a, 7r) = VAVijjV^VaVTi, (9) 
and where the source-dependent action Sj is 

Sj = 5gnjl + J^A^" + f]tlj + ^r] + J, a + J^vr + Sa. (10) 

The last term represents the ultraviolet regulating part of the action.^ The starting point 
to derive the Schwinger-Dyson equations (SDEs) is the following formal identity: 

j dfi (^A, ijj, Ip, a, tt) ^^^^^ exp {iSj} = 0, where (p = A,1p,^p, a,n, (11) 

from which the SDE for the propagators and vertices in momentum space can be obtained. 
For our purpose it is convenient to recall the SDEs for the scalar propagator and scalar 
vertex. The SDE equation for the scalar propagator is given by 

As\p) = -^ + Mp% (12) 
where the scalar vacuum polarization is 

Mp') = ^ f ^,Tr[S{k + p)Ts{k + p,k)S{k)] (13) 
Ja 

(see Fig. ^), S{k) is the full fermion propagator, and Ts{k + p,k) is the scalar fermion- 
antifermion vertex. The absence of kinetic terms for the a and n fields in the Lagrangian is 
refiected in the constant bare propagator — Go- 



^ In practice we use a regularization by means of a hard ultraviolet cutoff in the Euclidean 
momentum space integrals. Moreover the axial anomaly is neglected. 
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The scalar propagator is defined as the Fourier transform of the connected part of the 
correlator 

(0|T(a(x)a(y))|0)eonnected = / ^ e'^'^^^-^^^Aslfc) • (14) 

and the scalar vertex or Yukawa vertex Fs is defined as the Fourier transform of the ampu- 
tated vertex of the three-point correlator 

(o|T {ij,{x)MyMz)) |o)„ted = / e-^'=(^-^)+^^(^-^) 

X [zS{k){-z)Ts{k,p)zS{p)UzAs{k-p). (15) 
The full untruncated SDE for the scalar vertex in momentum space reads 

f d^r 

-iTsab{p + q,p) = {-i'i-)ab + J j^[iS{r + q){-i)Ts{r + q,r)iS{r)]a^ 

x{-iel)Kcd,ab{r,r + q,p + q) (16) 

(see Fig. |]), and the Bethe-Salpeter kernel K is defined as the two-fermion one-boson irre- 
ducible fermion-fermion scattering kernel. In the symmetric phase of the GNJL the pseu- 
doscalar and scalar propagators are degenerate, so are the pseudoscalar vertex and scalar 
vertex. 



III. SCALAR VERTEX IN QUENCHED LADDER APPROXIMATION 

In this section we discuss the SDKs for the scalar propagator and the scalar vertex in 
the well-known quenched-ladder approximation and introduce an approximation scheme for 
solving the SDE for the scalar vertex. The ladder approximation is obtained by replacing 
the Bethe-Salpeter kernel K by the one photon exchange graph. Furthermore the photon 
propagator is considered as quenched, i.e., vacuum polarization effects are turned off, and 
thus the gauge coupling does not run. In principle the Bethe-Salpeter kernel also contains 
scalar and pseudoscalar exchanges. One question is whether such exchanges can be neglected. 
It is beyond the scope of this paper to give a complete answer to such questions, since the 
answer not only depends on the short-distance behavior of the full scalar propagators and 
Yukawa vertex which is yet unsolved, but also on the representation of the chiral symmetry. 
In this respect, it is interesting to note that if one includes the ladder-like one-scalar and 
one-pseudoscalar exchanges in the truncation of the BS kernel K, and considers a chiral 
symmetry representation in which the number of scalars equals the number of pseudoscalars 
(thus both scalars and pseudoscalars in adjoint representation), then such contributions 
cancel each other exactly in the symmetric phase. Hence, provided scalars and pseudoscalars 
are considered both in the adjoint representation of the chiral symmetry, the neglect of scalar 
and pseudoscalar exchanges in the kernel K seems reasonable for the SDE for the Yukawa 
vertices Fg and Fp. 

The SDE equation for the scalar vertex in the ladder approximation can be written as 
rs(p + q,p) = l + tel [ 7^7^5(r + g)Fs(r + g, r)Sir)YDxAr - p) (17) 
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(see Fig. The SDE for the scalar propagator, Eq. ([T3|), is left unchanged. In the symmetric 
phase, the equation for the scalar vertex, Eq. (|T7|) , is a self-contained equation, if we note 
that in the Landau gauge the fermion propagator is S{p) = 1/p. The SDEs in ladder 
approximation for scalar propagator and vertex, Eqs. (|13D and (|T7D , have been studied 
extensively in the literature |7|,pJ],p!5|-p!7| , but mainly for the case of zero transfer boson 



momentum. 

In what follows, we present a method for solving the scalar vertex with nonzero boson 
momentum. The starting point is a general structure of the scalar vertex and pseudoscalar 
vertex. The scalar and pseudoscalar vertices in momentum space can be decomposed over 
four spinor structures with dimensionless scalar functions in the following way: 

Ts{k + q,k)= Ft\k + q,k) + (qk - kq) F^^\k + g, k) 

+ {k + q)Fi'\k + q,k) + kFt\k + q,k), (18) 



Tp{k + q, k) = (275) 



F^''\k + q,k) + [qk - kq) F!f\k + g, k) 



+ {k + q)Fl^\k + q,k) + kF^^\k + g, k) 



(19) 



where the scalar functions f/*^"* and F^^^ depend on the squares of the Minkowski momenta, 
{k + qY, k"^, q^. Thus by use of the notation F{k + q, k) we actually refer to a momentum de- 
pendence F{{k + qY, k"^, q^). Unless mentioned otherwise, our investigations will be focussed 
mainly on the symmetric phase of the GNJL model. In the symmetric phase, the equations 
for the scalar functions -F3 and F4 decouple from the equations for Fi and F2. Moreover, 
F3 and F4 do not contribute in scalar and pseudoscalar vacuum polarizations. In fact the 
integral equations for these functions are homogeneous ones and in the symmetric phase we 
can always take the solution F3 = F4 = which is a consistent one. Furthermore the scalar 
and pseudoscalar vertex functions coincide, i.e., f/^"* = F'f'\ i = 1,2. The consequence is 
that the scalar and pseudoscalar propagators are identical in the symmetric phase. So the 
problem is reduced to solving a coupled set of integral equations for two scalar functions Fi, 
F2, which has the form (after making a standard Wick rotation ]TB| ) 

F,{p + q,p) = 5a + \Y. j^^^^ j ■^K,,{p,q,r)F,{r + q,r), z = l,2, (20) 
where A = 3ao/47r, and 

^ / ^ 1 f^{p,q,r) . s 



( N 1 K{p,q,r){r^ + q-r) 



with 
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K[p, q,r)=pp-rq — p p ■ qq ■ r — 2p ■ qp ■ rq ■ r + 2p [q ■ r) 

+ 2{p ■ qYr"^ — 2p^g^r^ + p ■ rq^r'^ — p ■ qq ■ rr'^, (25) 

and / dQr denotes the usual angular part of the four- dimensional integration. 

The equations Eq. (^) are still very complicated due to the fact that the angular part 
of the integration cannot be performed in explicit form, since the angular dependence of 
the Yukawa vertex is unknown. Without any further approximations it seems impossible to 
solve the equations analytically. Our primary interest is the scalar propagator defined by 
the vacuum polarization Eq. (pIS]). The equation for the scalar vacuum polarization is 



A2 



\ r r (~\ \ } 

ns(g') = ^J J ^ [Mk, q)Fi{k + q,k) + A^^k, q)F,{k + q, fc)J , (26) 



where 

The method to tackle the angular dependence is to expand in terms of Chebyshev polyno- 
mials of the second kind Un{x), a method which was used before, for instance, in Refs. 
and PI. 

We define the following expansions: 

oo oo 

Fi{p + q,p) = ^ /„(p^g^)?7„(cosa), F2{p + q,p) = ^ g^)f/„(cosa), (28) 

n=0 n=0 

oo oo 

Ai{p,q) = J2(^n{p'^,q'^)Un{cosa), v42(p,g) = ^ 6„(p^g^)f/„(cosa), (29) 

n=0 n=0 

and for the kernels, Eqs. (0)-(|2^) 

oo 

K{p,q,r)= Knmiip'^,q'^,r'^)Unicosa)Umicosp)Ui{cos'j), (30) 

■n,m,l=0 

where 

p-q n P-r q-r 

cosa = , cosp = , cos 7= (31) 

pq pr qr 

(for the coefficients Knmi{p'^,q'^,f'^) see Appendix After that the angular integration 
can be done explicitly leading to an infinite chain of equations for harmonics /„(p^, q"^) and 
Qnip^iq^)- The important thing is that only the harmonics /o contains an inhomogeneous 
term in the equation for it (the constant 1 in Eq. (pOD), while other harmonics can be found 
iteratively once /o is computed. In other words, for the vertex function Fi the scale is set by 
the bare vertex, i.e., such a function has nonhomogeneous ultraviolet boundary conditions. 
For the vertex function F2 there is no such inhomogeneous term other than given indirectly 
by the coupling to vertex function Fi. 

We assume that the scalar vertex function Fi {p + q, p) depends only weakly on the angle 
between fermion and scalar-boson momentum p ■ q, so that an infinite set of equations for 
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/„ and gn is replaced by the equation for the zeroth-order Chebyshev coefficient function /o 
which we shall solve exactly. The main approximation is to replace the Yukawa vertex by 
the angular average of vertex function Fi, 

Ts{p + q,p) ^ 1 / ^F,{p + q,p) = l/ob',g'), (32) 

since Chebyshev polynomials of the second kind are precisely orthogonal with respect to 
such integration. Then we write 

/o(/, q') ^ F,^{p\ qyiq' - p") + F^^{p\ q^ip' - q'). (33) 

The functions Fir and F^y are respectively referred to as the IR channel (infrared), and the 
UV channel (ultraviolet). 

If the scalar vertex indeed weakly depends on angle between scalar-boson and fermion 
momentum p flowing through the Yukawa vertex, these channel functions should have the 
limits 

lim Ts{p + q,p) = 1 lim Fuv(p^ g^), (34) 
lim Ts{p + q,p) = 1 \im Fm{p'^,q'^), (35) 

i.e., the asymptotics of the scalar vertex are independent of the angle between p and q. 
Hence the UV channel contains a limit of the Yukawa vertex with the boson momentum q 
that is much less than both fermion momenta {q p), and the IR channel contains a limit 
of the vertex with the fermion momentum p that is much less than the boson momentum 
(g ^ p). The connection between the Yukawa vertex Fg and these two channel functions is 
illustrated in Fig. |. 

The expansion in Chebyshev polynomials is discussed in detail in Appendix Moreover, 
the error, due to our approximation, Eq. (0), in the computation of the scalar vacuum 
polarization, is estimated in that appendix. 

The zeroth-order Chebyshev or the two channel approximation of Eq. (|33|) gives the 
following equation for the vertex function /o(p^, q"^)'- 

/o(p^g') = 1 + a/ dr^No{r\p^)ao{r',q')fo{r\q^), (36) 



where 



Noir^p^)= ^ 2 + 2 (37) 



and 



/ 2 2\ -'- 

«o(r ) = 2 



^2-^)^(r^-g^) + ^%^-r^) 



(38) 



(see also Appendix 0). The equation for the scalar vacuum polarization in this approxima- 
tion Eq. (p6D takes the form 
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ns(g' 



47r2 



deao{k^,q')fo{k\q'). 



(39) 



With the Eqs. (0) and ( pSD for A^o and oq, respectively, and the definition of the channel 
functions Eq. (|33|), we get two coupled integral equations for the IR-UV channels: 



s t ^ 

Fir(s, t)^=^ 1 + xJdu -^^Fik{u, t) + xjdu ^Fm{u, t) + X J du 



2u - 1 
2m2 



FuvKt), (40) 



/u r 2,11 — t 
du —Aiu, t) + X du — Fuv(m, t) 
1st J 2su 



2u — t 

+ X I du Fuy{u,t), 



(41) 



where s = p^, t = q^, and u = r"^. For the vacuum polarization, Eq. (|5^), we obtain the 
equation 



ns(t) 



1 



47r2 



A2 



/It /" 2'U — t 

du—Fi^{u,t) + J du FuY{u,t) 



2u 



t = q\ 



(42) 



Using Eq. (^Tj) with s = provides a simple relation between the vacuum polarization and 
the UV-channel function 



ns(g' 



A2 1 



47r2 A 



i^uv(A',g'^ 



(43) 



which is different from the functional form proposed in Ref. 0. 

The integrals equations for Fir, and Fuv are equivalent to two second order differential 
equations with four appropriate boundary conditions. We get for the IR channel: 



ds2 

and for the UV channel 



s^^Finis, t) + 2s-^Fir(s, t) + ^sFir(s, t) = 



ds 



2t 



2s - 1 



s'^Fuv(s, t) + 2s— Fuv(s, t) + A^— Fuv(s, t) = 0. 
The infrared and ultraviolet boundary conditions (IRBC), respectively, (UVBC) are 



s'—Fi^is.t) 
ds 



0, 



s=0 



FuY + S — FuY 

ds 



1. 



(44) 



(45) 



(46) 



=A2 



Moreover we get a continuity and differentiability equation at s = t 



FiR(t,t) = Fuv(t,t), :r^iR(s,t) 

ds 



(47) 



s=t 
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The differential equations can be solved straightforwardly. The equation for Fir can be 
written as a Bessel equation, and the equation for Fuy as a modified Bessel equation. The 
general solutions of the differential equations are 



^2Xt/s 



(48) 
(49) 



with Ji and Yi the Bessel functions of first and second kind, respectively, and where are 
modified Bessel functions and u = \/l — 4A, see also Eq. (|). We note that since Eqs. (^4|), 
(^) are scale invariant their solutions are functions of the ratio s/t and the scale invariance 
is violated by the UV boundary condition ( ^61) only. 

The IRBC for the IR channel requires C4,{t / A"^ , uj) = 0, since the Bessel function Yn{z) is 
irregular at 2; = 0, the other coefficients are fixed by the remaining three boundary conditions 
and the solutions are 



c,{t/A\u) = f^Z-\t/A\u), 
2 smuJTT 

C2(t/A^cu) = cl(t/A^-cu), 

C3(^/A^a;) = Z-l(t/A^cu), 



where 



and 



Z{t/A^,uj) 



TC 



2 sincjvr 



7(cu)G(t/A2, _ ^(_^)G'(t/A2, cu) 



(50) 

(51) 
(52) 



(53) 



7(^) 
G{t/ A\u}) 



Ji(V2A)/:(V2A) + J((V2A)/^(V2A) 



1 U 

2 V A2 



2At/A2 - J2\t/AH'^ J2\t/A 



Summarizing, the solution of scalar vertex in terms of channel functions is 



,2\ 1/2 



A2 



TT 



X 



2 sincjTT 



-z 



-1 / ? 



Jl 



l2Xp^ 



VA2' 

/2V 

p2 



UJ 



2\ 1/2 



-7(-^)/. 



/2Ag2 



It is easy to verify that at zero boson momentum one gets 



rs(p,p)=Fuv(p^g' = 0) 



2 P_ 

1 + uj \A^ 



2 X -il-u;)/2 



(54) 
(55) 



(56) 



(57) 



(58) 
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which coincides with the zero transfer vertex of Refs. |]7|JTl[] . In the hmit of pure N JL model, 
tto — >■ (u; ^ 1), the vertex is equal, of course, to the bare vertex, = 1 = Fuv = -Fir. 
To study the vertex at critical gauge coupling = etc {u = 0) we expand the Bessel 
functions in small uj using the following property of the modified Bessel functions: 



^ /o(x(0)) - ujKo{x{Q)) + 0{u:^ 



where x{uj) (x yl^^uj^. Then the expressions for Fir, Fuv, Eqs. 
form at ao — * etc: 



where 



2 \ 1/2 



A2 



,0 



1/2 



p2 



2q' 



2 \ 1/2 



- £2^0 



(59) 

take the following 



(60) 
(61) 



[^i^oix) + eixKi{x) - 62/0(3;) - e2xlQ{x)] 



X 



g2/2A2, 



(62) 



61 = Vl/2(Ji(Vl/2)/o(\/V2) + J((Vl/2)/o(Vl/2)), 

62 = /l72(Ji(/l72)ir;(/l72) + J[i^)K,{^)), 



(63) 
(64) 



and Ki is the modified Bessel function of the third kind. Note that when we expand the UV 
channel function, (p^ ^ g2-j^ (|6l|) , we get 



^uv(p',g') =2 



A2 



-1/2 



£3 - 2 + ln(p2/g2^ 
es - In(g2/A2 



(65) 



where 



€3 = 2(1 - 7) + 3 In 2 - 2- ^ 3.2 

ei 



(66) 



and 7 is the Euler gamma. The expression (KH) is of the same form as obtained in Ref. 
(their formula (2.88)). 



IV. THE SCALAR PROPAGATOR 



In the previous section we obtained an analytical expression for the scalar vertex by 
assuming that the vertex only weakly depends on the angle between boson momentum and 
fermion momentum. The expression for the scalar vacuum polarization in such an approxi- 
mation takes the form of Eq. (|42|) . The main object of investigation is scalar compositeness 
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near the critical line. In the neighborhood of the critical line, the tendency of fermion- 
antifermion pairs to form bound states under the influence of strong attractive four-fermion 
forces becomes apparent. 

Since we are in the symmetric phase of the GNJL model (there is no dynamical mass), 
the only important variable is the scalar boson momentum over cutoff t = g^/A^. Equa- 
tion (^21) shows that the vacuum polarization depends only on the UV channel with fermion 
momentum at A^, Eq. (^Hj). Substituting the expressions obtained for the vertex function 
Eq. (^9]) in the equation for the vacuum polarization Eq. (53), we obtain 



A2 1 



47r2 A 



7(cu)i/(gVA^ -uj) - ^{-uj)H{q^lK\ 00) 
, 7(c^)G(gVA2, -co) - ^{-uj)G{q^ / A\ co) 



where 7 and G are given by Eqs. 



and 



and where 



t 



Hit/A\uo)^-J- 



(67) 



(68) 



This expression is valid along the entire critical curve in the symmetric phase. Note that 
the vacuum polarization is symmetric in lo which means that we can analytically continuate 
it to values ao > etc- 

Equation (^) is a rather complicated expression, and we will first investigate some 
specific limits: 

(A) The pure NJL limit, where the gauge interaction is turned off, i.e., the case where 
ao = 0, thus a; = 1. 

(B) The behavior at the critical gauge coupling ao = Oc, thus lo = 0. 



(C) The behavior of lis for c^o > etc? ^ 
across the critical curve at ao = Oic 



1, i.e., analytic continuation 



(D) Asymptotic behavior of Ilg, the infrared behavior g^/A^ <^ 1 in such a manner, so 
that {q^lA^Y > gVA2. 



A. The pure NJL limit 



The pure Nambu-Jona-Lasinio limit is the case where the gauge interaction is completely 
turned off, possible bound states in such a model are purely due to the four-fermion interac- 
tion. So tto 0, and therefore u; — >• 1. The scalar vertex at cto = is just equal to the bare 
vertex, Fs = 1 = -Fuv = -Fir. The pure NJL limit can be correctly obtained from Eq. ( |67D 
by making expansions of the Bessel functions keeping sufficient number of terms, and then 
performing an expansion in (1 — lo) to zeroth order. Hence with ao = 0, = 1 in Eq. (p7D, 
we get 



ns(g^ 



A2 
47r2 



'+2A^^"U -4A^ 



(69) 



The log term is a consequence of the "hard" fermion loop which also ruins the renormaliz- 
ability of the pure NJL model. 



12 



B. The a boson at critical gauge coupling 



At critical gauge coupling ao = ac {uj = 0) is the onset of the scalar compositeness 
origining purely from electromagnetic forces. The expression for Fuy, Eq. (0), has the 
form at ao = 



P 



2\ -1/2 



e2/o(y) - ei^o(y) 



e2/o(x) - e2xli{x) - eiKolx) - eixKi{x) 



(70) 



where x 



j2h? and y = \Jq^/2p^. This gives 

r. S[e2xli{x) + eixKi{x)] 



4 + 



€2/0(3;) - e2xli{x) - eiKo{x) - eixKi{x) 



(71) 



C. Analytic continuation across the critical curve 



Since the expression for the scalar vacuum polarization is symmetric in u, Eq. ([67|), it 
can be analytically continuated to the values ao > c^c- This holds in replacing u hj iu in 
Eq. (|67D, where 



= Jao/ac- 1. 



(72) 



D. Asymptotic behavior of the a boson vacuum polarization 

The asymptotic behavior < < 1 so that {q^ / K^Y ^ g'^/A^ can be obtained by 
considering first the <^ A^ limit of Z, Eq. (|53|): 



7r 



2 \ 1/2 



where 



2 /i(-c<;)(l -cu) 



uo In a/2A, /i(c<j) = 



7(^)2" 

rri -cu) 



(73) 



(74) 



Then the UV-channel function with fermion momentum = can be expressed, in this 
limit, as 



Hence 



ns(g' 



2 2cj , , , cj, /A2\ ^ 

^ + ^(l-coth,), ,= -ln - +^(.). 



A^ 

477^ 



g,{L0) (1-^2)2 



'1 — cothy) 



q^ « A^. 



(75) 



(76) 
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However, the above expression does not reproduce the correct leading term of the NJL hmit 
(cu — i> 1), for that we should use the expression (|67D . In order to obtain such an expression 
which contains properly the pure NJL limit, we have to keep more terms in the expansion 
of Z. We then get for the scalar vacuum polarization 



ns(g') 



A2 
47r2 




B{u) ) + A{uj)^- + O ((gVA^)^i + O ((gVA^) 



A2 



where 



A{uj) 



2g,{uj){l-uy 



B{uj) 



IQoj 7(-c<;) r(2 - cu) f\" 
^ 7(cj) V{2 + uj) I 2 



(77) 



(78) 



For Eq. (|77D to be valid, it is assumed that > 1/2. It is straightforward to check that this 
expression satisfies the NJL limit {uo 1). 

This is a suitable point to refer to Appendix ^ for an analysis of the reliability of Eq. (ffTD, 
where it is shown that the leading and next-to-leading terms in g^/A^ of lis (i-6., the first two 
terms on the right-hand side of Eq. (^)) are indeed correctly obtained by our approximation, 
see Eq. ( |A35| ). 

The four specific limits of the scalar vacuum polarization described above are very useful 



for studying the resonance structure of the bound states which will be done in Sec. 0, and 
for the study of the CPT which will be done in Sec. |Vli| . But first we shall compare our 
result for As with that one obtained by other authors in the next section. 



V. COMPARISON WITH EARLIER WORKS 



In this section we discuss the earlier work of Appelquist, Terning, and Wijewardhana 
||Tl| and related work based on that by Kondo, Tanabashi, and Yamawaki on the scalar 
composites in the GNJL model. The method used by these authors to solve the coupled set 
of scalar vertex and scalar vacuum polarization is the following. They consider a Taylor- 
series expansion about g = of the scalar vacuum polarization ns(g^). In the ladder 
approximation, such a series has the property that the nth derivative can be written as 



n— 1 



I A {2tt) 



m=0 



' Qm \ Y Qn-m ^ 

^rs(fc + q,k) \ ^rs(A;, k + q)- 



■ (79) 



Their basic assumption is then that derivatives of the scalar vertex Fs can be neglected with 
respect to Fs for small ao. 



9"Fs(A; + g,A;) 



dq^- 



5=0 



1 «o 



Ts{k,k). 



(80) 



Subsequently the Taylor series is resummed using the assumption stated above to obtain 



ns(g^ 



Tr 



/a (27r)4 

which yields in terms of Euclidean momentum. 



Fs(A;,fc)iFs(fc,A;)J- 
K k — q 



(81) 
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where 



1 



4vr2 



1 ( q'^Y q^' 



a 



1 



9c 



UJ 



(82) 



(83) 



How does their result compare to ours? From our expression for the asymptotic behavior 
of lis, Eq- (|T7|), and the result obtained in Refs. |Tl| and Eq. (p2D, the leading power 
of momentum is the same, namely, (g^/A^)'^. However, the tu-dependent factors in front of 
the leading and next-to- leading powers are different. At the same time the pure NJL limit 
is obtained correctly in Refs. |]lT| and 0. The differences are rather small for values of lo 
close to 1, i.e. the coefficients are 



A{uj) = a{u;), B{lj) = b{u;) 1 + 0{{1 - uf) 



(84) 



Hence the approach to the NJL point u = 1 oi both approximations is equal. However, for 
smaller values of uj the coefficient B obtained in the present paper starts to deviate from 



the one obtained in Ref. |Tl|. Then for such values of b{u) > B{u). 

What is the origin of this difference? The first point is that the expression derived in 
Refs. |]TT|] and [0 is valid for ao not too large. Secondly, from their answer Eq. it is 
clear that a Taylor series of the scalar propagator about g = is not well defined due to the 
noninteger power behavior for values < u < 1. This is reflected also in their assumption 
regarding the derivatives of the scalar vertex at g = 0, Eq. (pOl). The expression for the 
scalar vertex obtained in Sec. |TTl| shows that in general, for < u < 1, the assumption (|80D 
is not true. Such derivatives of Fs are singular at g = due to the fact that they depend on 
noninteger powers of g, which can be seen from Eq. 

The scalar vertex for small g <^ p is of the form 



-1/2+1^/2 



q^ 



P_ 

A^y \{i+u)^ p 

2\"^(-o;)F(l-o;) 



uj) 




7(cu) F(2 + a;) 



4 4(l + cj)A2 

, /.2 





+ O {{qVp^riqVA^r) + O ((gVp')'+" 



+ 



^5) 



which is consistent with Eq. (^) for p"^ = A^, because of Eq. (|^). Of course at g = our 
result coincides with that of the other authors. But for nonzero g this expression clearly 
shows that for < u < 1 the vertex contains noninteger powers of g. Hence the assumption 
made in Refs. |]7|,p!T|| is not true in general, since higher derivatives of the scalar vertex with 
respect to g are singular at g = 0. 

In what follows, we discuss the results from the viewpoint of the renormalizability of the 
GNJL model. The renormalization is performed by a suitable redefinition of the composite 
or auxilary fields a and tt: 



0"R 



-1 



IT. 



(86) 
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The renormalization factors of the scalar and pseudoscalar fields Zjf\ respectively, Z^'^^ 
can be chosen to coincide both in the symmetric and broken phase, and the renormalized 
auxilary fields an and tt^j define the renormalized scalar propagator 

-2 



As(g) 



and the renormalized scalar vertex 

rf(p + g,p) = 4^)rs(p + g,p). 



(87) 



(88) 



In order to renormalize the scalar propagator and Yukawa vertex simultaneously (see also 
Eqs. ( |77D and (p5|)), the wave function renormalization factor at some arbitrary renormal- 
ization scale fi should be of the form 



/ 2 \ (l-'^)/2 



(89) 



Freedom in the choice of renormalization scheme allows us to take the factor Z defined in 
Eq. ( p3D as the wave function renormalization factor, since 



Z('^) = Z(/iVA^a;) 



vr 



7(^) (1 + ^) 



2 sincjvr 2 Ff 1 — u) 



-uj/2 



A2 



(l-cv)/2 



(90) 



Hence it follows that four-fermion scattering amplitudes, for instance one-scalar exchange 
amplitudes are renormalization group (RG) invariant, i.e. 



(pi + q,pi)A^{q)T^{p2,P2 + g) ~ Fs(pi + g,pi)As(g)Fs(p2,P2 + q)- 



.Ri 



(91) 



Consider the case where pf, P2 ^ Q'^, so that the scalar vertices are described by the UV 
channels, and suppose that we are sufficiently close to the critical line 



Ago = 9c- go ^ ij^ 



< g', 



then from Eqs. (p 



and ([77| ) the scalar propagator has the asymptotic behavior 

47r^ 1 
' A2 B{uj) \ q^ 



As(g) 



(92) 



(93) 



Such a specific power-law behavior for the scalar propagator is essential for the renormaliz- 
ability of the GNJL model as is shown in Refs. 0-0] • 
Thus, from Eq. (RUf) and Eq. (B^), we get 



Af (g) oc - 



UV 



(p^g^) 



1/^ 

2 / 



oc 



p" 



7(-cu) F(l - cj) /V 



(94) 
(95) 



With these expressions, it is straightforward to check that Eq. (^Tj) is indeed independent of 
A and ^. Hence the renormalization of the auxilary fields a and tt, Eq. (0), simultaneously 
renormalizes the Yukawa vertex and the scalar propagator. 
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VI. LIGHT SCALAR RESONANCES NEAR CRITICALITY 



In this section we discuss the behavior of the scalar propagator near the critical line in the 
symmetric phase 5^0 < Qc- In the symmetric phase the scalar and pseudoscalar composites, 
the a and tt bosons are degenerate. Near the critical curve, a combination of strong four- 
fermion coupling and gauge coupling will tend to bind fermions and antifermions into these 
scalar composites. Since the chiral symmetry is unbroken the a and vr bosons decay to 
massless fermions and antifermions. Hence the scalar composites are resonances which are 
described by a complex pole in their propagators. The complex pole determines the mass and 
the width of the resonances. In what follows, we redo the computation of the complex poles 
of the a boson which was performed by Appelquist et al. in Ref. |Tl[ using the expression 



for As, Eq. (|67|) , obtained with the two-channel approximation of the Yukawa vertex. The 
expressions obtained in Sec. |rV| for ns(p^) in various regimes are rotated back to Minkowski 



momentum p\j exp^—iir) . Then the complex poles are given by 



Pm =plexp{-ii 



^s\P 



A' 



M 



(96) 



We can also parametrize the location of a pole by a mass and a width, i.e., p^exp {—i6 
[M^ - {i/2)T^f, which yields 



Po 



1 + cos 6^ 



nl/2 



2 sin 6* 
1 + cos 6* ■ 



(97) 



If 9 is small, then T^/M„ ^ 9. 

Near the Nambu-Jona-Lasinio point (oq = 0) our expression for the vacuum polariza- 
tion coincides with that obtained by Appelquist et al., Eq. (0), and we get the following 
equations for the resonances: 



cA2 _ ln(AVp2) + 3/2 _ 2(1 - ^p) 

, c 



Po 



cos 9 



9o 



(98) 



and 



n + 9 



H^Vp'o) + 



tan^. 



(99) 



If now qq is tuned close enough to the critical value Qc = 1, so that In 1/c ^ 1, the solution 
is approximately 



pI 



2(1-90) 



and we find a narrow width 



9 



go\n [5-0/2(1 - go)] 



TT 



ln[V2(l-5o)]' 



A^ 



(100) 



(101) 



These results are nothing else than the familiar NJL results. For intermediate values of the 
gauge coupling, < ao < etc (0 < < 1), we assume the poles of As are small, Po/A <^ 1, 
so that 
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(Po/A)" « 1. 



(102) 



Then, from Eq. (|77|) we get the following equation for the real part of the pole: 

: 

The equation for the imaginary part reads 



9o 9c 







sma; 



{0 + n), 



(103) 



(104) 



where B{uj) is given by Eq. 
The solution is 



7r{n — uj) 



(105) 



and n is odd integer, so that cosco'(6' + vr) = —1, thus 



Po ~ A 



goBiuj) 



l/2u) 



(106) 



Hence 9 is only small if ti; ~ 1 for n = 1. The result obtained in Ref. [|lT[], see Eq. (|83|) , gives 
a mass 



Po ~ A 



'1 - 



l/2ij 



1 



(107) 



The pole obtained in Ref. |Tl|] is of the same order as Eq. (|106|) for values of u close to 1 
(see also discussion in the previous section). For more intermediate values of u the poles 
obtained in our approximation are somewhat bigger, since b{uj) > B{uj) for < < 1. The 



quantitative difference between the result of Ref. |TI| and that obtained in this paper for 
resonance structures are visualized in Fig. In Fig. ^ the imaginary part of As given by 
Eq. ( p7[ ) and Eq. (|T^) is plotted versus p/M where the tuning of the four-fermion to the 
critical line is go/ Qc = 0.999, and M/A = (1 — g^/ Qc)^/"^. From Fig. |^ it is clear that the 
position of the peak of the resonant curve is slightly shifted to the right in our case at a 
fixed ratio ao/«c ~ C^(l) (intermediate or small values of a;), while the width over mass 
ratio remains comparable. Near the pure NJL point ~ 1 both results coincide. 

As was pointed out in the previous section, and following from the restriction Eq. ( |102D , 
these results are only valid for a^/ac small. For larger values of the gauge coupling, u — >• +0, 
the widths become larger, and Eq. (|102[) is no longer satisfied. This can also be seen in Fig. 0. 



As the ratio ao/ac increases the width increases too, and the position of the peak becomes 
more difficult to define. 



VII. THE CONFORMAL PHASE TRANSITION (CPT) 

In this section we analyze the scalar composites near the critical gauge coupling = 
with the purpose of investigating the conformal phase transition. The conception of the 
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CPT was introduced and elaborated recently in Ref. [0]. It embodies the classification of 
specific types of phase transitions. The main feature of the CPT is an abrupt change of 
the spectrum of hght excitations (composites) as the critical point is crossed, though the 
phase transition itself is continuous. This is connected with the nonperturbative breakdown 
of the conformal symmetry by marginal operators {{ipip)'^ + {ipi'^^ipY ^^e model under 
consideration), which was illustrated in Ref. [|14| by a study of the effective potentials in 
Gross-Neveu and GNJL models and quenched QED4. 

In the previous section we encountered a no-CPT, cr-model-like phase transition for values 
of ao < Uc [0 • The masses of light excitations^ are continuous functions across the critical 
curve; there is no abrupt change in the spectrum of light excitations. In the broken phase 
the TT boson becomes a massless Nambu-Goldstone boson, while the fermion and a boson 
acquire a dynamical mass which is small with respect to the cutoff A near criticality. 

The pole at the critical gauge coupling = Uc is determined by Eq. (^), which is rather 
complicated but we only need to study the IR limit, so we assume that the pole is small 
Pq <^ A^. The infrared limit obtained from Eq. (fflD is 



ns(g' 



A2_ 

47r2 



4 + 



16 



In(gVA') - ^3 



+ 0(gVA2ln(gVA^ 



(108) 



We then find zeros of Ag^{p) at 




16(ln(AVpg) + es) 



(109) 



0^9 + 71. (110) 

Thus 6 ~ — vr, and from Eq. ( |109| ) it is clear that ii go < Qc = 1/4 both terms on the 
right-hand side are positive, and there is no solution for the pole with po/A <^ 1. Hence if 
there is a pole it will be heavy, i.e., po/A = 0{1). Therefore at ao = etc no light resonances 
are present in the spectrum for 5^0 < fl'c = 1/4. The imaginary phase 6 approaches — tt which 
means the heavy pole occurs at "Euclidean" momentum, a sign of tachyonic states. 

The statement above can be made more explicit. If we analytically continue the scalar 
propagator to the values of ao > ac, then we end up in the "wrong vacuum" and we should 
get tachynonic states. In the broken phase (ao > etc), a chiral symmetric solution still exists, 
but it is unstable. The tt and a bosons are tachyons for such a solution. The unstable 
symmetric solution is obtained by analytic continuation of the solution in the symmetric 
phase across the critical curve (at ac)- The scaling law is determined by the UV properties 
of the theory and therefore the scaling law of the tachyonic masses is the same as that of 
the fermion and a-boson mass in the broken phase. 

Tachyons are described by imaginary mass < 0. This means the scalar propagator 
must have a real pole for Euclidean momentum. If the pole po is small, po <^ A, we 
analytically continue Eq. (^) to ao > «c by replacing u by iu, 



^The scalar composites are light resonances. 
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u ^iu = iV4:X - 1. (Ill) 

We then obtain 

= ^ ' y=2'"[-^)^-'p(-''>' (112) 

where 

0(.^) ^ ^ In ^^(^'^\^; + ^'^\ - In V2X. (113) 
The tachyonic pole is then given by the zero of the equation 

A2 



which gives 



+ ns(p^) = o, (114) 



— = expl — — -- + 0(z/)l, (115) 

where n is a positive integer and 

P = tan-i -. (116) 

The tachyon with largest po in the physical region < K corresponds to n = 1. If we now 
consider the limit i/ — > (A —> 1/4), we get 

'^^90 ,2^ _ 1 1, o ^2 , ,^,.,2x 



In this case, 



(3^ 0(i.^)^l-7--ln2-^ + O(z/^). (117) 

9o - 1/4 2 ei 



which is proportional to the well- know scaling law of quenched QED. Thus the scalar propa- 
gator giving the tachyon pole equation, Eq. (|118|) , reproduces the scaling law with essential 
singularity, which is another confirmation of the CPT. 



VIII. CONCLUSION 

In this paper we studied the scalar composites near criticality in the GNJL model. We 
obtained an analytic expression for the scalar propagator describing the composite states 
which is valid along the entire critical curve of the GNJL model. We presented a method for 
solving the Yukawa vertex in the GNJL in the quenched-ladder approximation. The crucial 
assumption was that such a vertex depends only weakly on the angle between a-boson 
momentum and fermion momentum. The method presented here incorporated the infrared 
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boundary condition in a more natural way than previous attempts in this direction. Also the 
observation that derivatives of the Yukawa vertex are singular at zero a-boson momentum 
transfer is a warning that derivative expansions and Taylor series could fail. Moreover this 
is reflected by the property of the scalar composites having noninteger power-law behavior, 
which means that, although these states are tightly bound, they are not pointlike. 

The conclusion of the comparison of the method presented here and work done previously 
on the (T-boson propagator is the following. Qualitatively the results obtained by Appelquist 
et al., Kondo et al. and by our approximation are in agreement. Both methods yield a renor- 
malizable a-boson propagator and Yukawa vertex near criticality and find light resonances 
for gauge coupling < ao < Oc- Quantitatively, at a fixed value of the gauge coupling, the 
scalar composites computed in the our case are slightly heavier with comparable width, as 
is illustrated in Fig. |^. 

In addition, the scalar composites propagator were examined for values of the coupling 
near the critical gauge coupling ac- Near the critical line = ac, go < \ the conformal phase 
transition is encountered and the spectrum of light excitations (resonances) in the symmetric 
phase disappears. Moreover the well-known scaling law with essential singularity, which is 
characteristic for the CPT, was recovered by analytic continuation of the a-boson propagator 
across the critical curve at ao = ac- 
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APPENDIX A: ANALYSIS OF THE CHEBYSHEV EXPANSION 



In this appendix we discuss the validity of the zeroth-order Chebyshev expansion for the 
Yukawa vertex function Fi introduced in Sec. |T|. The problem of angular dependence in the 
SDKs for the Yukawa vertex functions Fi and F2 is replaced by an infinite set of Chebyshev 
harmonics. Subsequently this set is truncated to the lowest order harmonic of Fi, which is 
the only harmonic having nonhomogeneous ultraviolet boundary conditions because of the 
presence of the (angular independent) inhomogeneous term 1. 

As mentioned previously, the method of using expansions in terms of Chebyshev polyno- 
mials Un{x) (of the second kind) was used before [^|19|. These polynomials are orthogonal 
with respect to the angular integration / dQ. In the analysis of BSEs in Ref. ^ a CP in- 
variant Chebyshev expansion was used, which has the nice property of keeping only even 
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terms in the expansion. However, we use a slightly different expansion (not explicitly CP 
invariant) which has the disadvantage of also including odd terms in the Chebyshev expan- 
sion, but the advantages that the integral equation for the zeroth order harmonic is more 
friendly and the zeroth order harmonic coincides with both the large fermion momentum 
limit (p^ ^ g^) as well the large boson-momentum limit (g^ ^ p^) of the Yukawa vertex, 
see Fig. ^. 

Thus the vertex functions satisfying the SDEs (^) are expanded in the angle between 
fermion momentum p and scalar boson q, i.e., p ■ g, in the following way: 



Fi{p + q,p) = J2 fn{p'^,q^)Un{cosa), 
1 



n=0 

oo 



n.=0 



(r — pY 



J2N^{r',p')Un{cos(3), 



n=0 
oo 



n=0 



F2{p + q,p) = ^gn{p^,q^)Un{cosa), (Al) 

(A2) 
(A3) 



where 



cos a = , cos p 

pq 



A2{r,q) = ^ 6„(r2,g^)f/„(cos7) 

n=0 



p ■ r q ■ r 

cos 7 



pr 



qr 



(A4) 



The vertex functions and kernels Ai and A2 were defined in Eq. ([T8|), respectively, Eq. 

g 



The coefficients A^„, a„, and 6„ are 



Nn{r\p') 



and 



ao(r^g2) 
an(^^^g^) 



2-^-]9ir'-q') + ^9{q'-r' 



(r'^-q^) 



Oir"^ — q^) f q\^ 9{q^ — r"^) /r' 



n > 1, 



(A5) 

(A6) 
(A7) 



and 



6o(r^g2) 
h{r\q^) 



^2 ^V^ -9 ) + 



-9{r' - r) + 



r^(r^ — 



eiq^ - r^) 



q6 



2 2\ 

— r ) 



2\2 



hn{r\q') = {-!)' 



(^2 _ ^2) 



Oir"^ — q^) /g\" 9{q^ — r"^) f 



r , ^ 



(A8) 
(A9) 

n > 2. (AlO) 



The equations for the scalar vertex functions ( PP] ) and scalar vacuum polarization 
are expressed in terms of an infinite set of equations between the harmonics. Hence 

ns(g') = j^JdeY: h(A:^ q')fn{k', q') + 6.(A;^ g2)^„(A;^ q')] , (All) 

n ^1=0 
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and with Eqs. (pT]) and (p^) using Eqs. and ([A3|) , we get for the harmonics of Fi 

A 



^ ' 



oo oo 



X E Y.^irnn[ara{r'^,q^)fn{r'^,q^) + hm{r'^,q^)gn{r'^,q^) 

m=0 71=0 



where 



(A12) 



Cimn = - / d7 sin^7?7;(cos7)f/m(cos7)f/n(cos7). 
vr J 



(A13) 



In the derivation of Eq. (|A12| ) use has been made of the fact that 

1 cos(a-7) 



1 



7rsm7 



dcosa 



dcos/9?7m(cosa)?7„(cos/3) = 5m,n 



cos(a+7) 



[4(cos7) 
n + 1 



(AM) 



The symmetric index Cimn can be calculated using product properties of the Chebyshev 
polynomials, giving 



Co 



I 

I 

m—n\ ^2k-\-l,n-\-m-\~2 

/ = 0,1,2, 

k=0 



(A15) 
(A16) 



Then the first two equations for the coefficients /n(p , q ) of vertex function Fi read 

\ OO 

/o(s, t) = 1 + X du No{u, s) [«m(M, i!:)/m(M, i^) + bm{u, t)g^{u, t)] , (A17) 



m=0 



X r °° 
/i(s,t) = - duN,{u,s)Y: 



am+l{u,t)fm{u,t) + am{u,t)fm+l{u,t) 



+ bm+l{u,t)gm{u,t) + bm{u,t)gm+l{u,t) 



where we have introduced the variables 



2 J 2 2 

s = p , t = q , u = r . 



(A18) 



(A19) 



In principle there is an equivalent set of equations for the coefficients gn of F2. However 
we did not succeed in finding explicit expression for these, due to our inability to compute 
explicitly the angular integrals given by kernels 7^21 and K22 of Eqs. ( P5D and(^). The 
problem is to compute the integrals 



dVl„ f dflr 



J ^U4cosa)K2,ip,q,r), J^J ^[/„(cos«)K2i(p, g, r). (A20) 
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The main approximation used in this paper is Eq. (p2D, the replacement of the Yukawa 
vertex by the zeroth-order harmonic of Fi. In what follows we estimate the error made by 
such an approximation. We define the error E{q^) in computation of the scalar vacuum 
polarization Eq. ( |A11|) as follows: 



ns(g' 



dk'a,{k\q')U{k\q') + E{q') 







.n=l 



n=0 



(A21) 
(A22) 



For an estimation of E{q^) we need to know more about the harmonics /„, n > 1, and 
> 0. The solution to these harmonics is assumed to be governed by the harmonic 
/o only, since the integral over the harmonic /o acts as the largest inhomogeneous term in 
the integral equations for the higher order harmonics. So the equations for the higher order 
harmonics are approximated by 



gn{p\q') ^^J<^^" J^J ^U^i^osa)K,,ip,q,r)Ur\q'), 



X 



(n + 1) 



(A23) 
(A24) 



Unfortunately there is no explicit expression for Eq. ( |A23| ) for the reason described above. 
However it is possible to approximate the angular average by considering either one of the 
three momenta in K2i{p, q, r) to be much smaller than the other two. Then the dependence 
on one of the three angles between the momenta is lost, and the integration can be performed 
explicitly. The result for the lowest harmonic of F2, i.e., given by Eq. ( |A23| ) , is 



go{s,t) 

9ois,t) 



— / ' 
12 Jo 

X ft 





[u,t) + 






-/ 

12 Js 




>,t) + 






-1 

12 Jt 



12 Jq 

and for Eq. (|A2l ) with Eqs (^) and (^) 



^ . u-t V? ^ , , X f*-^ u-t fs ^ , 
X u-t [si 



4 



— —\\—Fxjy{u,t) 
V M 



, X f^^ U-t lu^^ , , X u-t , . 

/i(s,t) ^ -- / du — —\ —FiK{u,t)-- du W-Fuv(M,t) 

4 JO st \ st 4 Jt us \ s 



X u-t [si 

du 



9 \/ 9 
V M 



Fuv{u,t), 



(A27) 



(A28) 
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where s = p , t = q , and we have used Eq. (|33|). These equations can be analyzed in detail 



once the solutions for the channel functions, Eqs. ( ^6|) and (57), are known. But for obtaining 
the asymptotic behavior of the harmonics Qq and /i it is sufficient to use the asymptotics of 
the channels, i.e., take Fi^ip'^, q^) Fir(0, q^) oc (gVA2)^/2-V2^ Fuvb^ q^) ^ i^uv(p^ 0) cx 
{j? I t^Yl"^^^!"^ . This gives for the harmonics 

g^{v\ q') oc x\fMp^ 0), q') (X X^F^p', 0), / » g^ (A29) 

go{p\q^) oc A^FiR(0,g2), h{p\q^) oc A^FiR(0,g2), ^^gg^ 



The above equations give the leading behavior (in either q/p ^ 1 or p/g ^ 1) of the har- 
monics ^^o; /i in terms of /o up to some A-dependent factor, which is 0{1) (thus nonsingular 
in A). Furthermore, from Eq. ([A24|) we get the relation 



fn+i{p,q^) ^ _q 
fn{p'^,q'^) p' 



p" > q\ 



/«+i(p ,g ) ^ _P 
fn{p^,q^) q'' 



q^ > p\ 



(A31) 



and we assume a similar relation to hold between the harmonics Qn+i and Qn- Thus the 
series 



S{p^,q^)= [(^m{p^,q'^)fm{p'^,q^) + bm{p'^,q^)gm{p'^,q^ 



m=0 



(A32) 



which occurs both in the equation for /q, Eq. ( A17 ), and for lis, Eq. ( All ), will be rapidly 
converging for either p"^ ^ or p^ <^ q^, since 



Qn+l(p^g^)/n+l(p^g^) _ min(p^g^ ^ 

«n (P^ ,q'^)fn{P^,q'^) max (p2 , g2 



«1, 



(A33) 



and again a similar equation for the part containing the harmonics Qn- At p"^ = q^ only three 
terms of the series S, Eq. (|A32|) , contribute, since anip^.p^) = for n > 1 and 6„(p^,p^) = 
for n>2. Hence, a straightforward approximation for the series S is 



S{p\ q') ^ ao(p^ q')fo{p', q') + O {a^{p\ q')f,{p^ q')) + O {bo{p\ q')g,{p\ q' 



(A34) 



supporting Eq. (|32D . With the expressions obtained for /i and go, Eqs. ( |A29|) and (|A30D , 
the leading term of the error E defined in Eq. ( |A22| ) can be estimated. The leading term of 
the error E[q'^) is given by 



E{q' 



47r2 



a^{e,q')Me,q') + bo{e,q')go{e,q'' 



A2 



A J dk'^ ^Fuv(A:^, g^ = 0) + next-to-leading 



A' 



;i + o(A)) I, 



2 \ cj/2+1/2 



:i + c^(A))^ 



(A35) 
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where we have kept only leading terms and for Fuv(p^,0) = Fi{p,p) given by Eq. 
Recall that uj = vT^-4X. The estimation of Eq. ( A35| ) can be checked more explicitly by 



using the solutions obtained in Sec. fT| for the channel functions Fir, Fuv, Eqs. (|^) and 

Eq. ( |A35| ) shows that when X = 0, uj = 1, the error E vanishes, and when uj < 1 clearly 
the terms in the error can be neglected with respect to the first two terms on the right-hand 
side of Usiq'^), see Eq. ([77[). Thus this analysis supports the assumption Eq. ( [5^ ) made in 
Sec. m and the error E contributes only to next-to-next-to- leading order in g^/A^. And 
therefore we may conclude that our approximation gives correct leading and next-to-leading 
behavior of ng(g^). 
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FIGURES 




FIG. 2. The SDE for the scalar propagator As(p). 
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FIG. 3. The SDE for the scalar vertex Ts{p + q,p)- The shaded circle with the 2 represents the 
two-fermion one-boson irreducible fermion-fermion scattering kernel. 
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FIG. 4. The SDE for the scalar vertex in the quenched-ladder approximation. 
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FIG. 5. The two-channel approximation. 
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golQc = 0.999 




p/M 

FIG. 6. A comparison of our result with Appelquist et al. Graphs of Im(As(p)) for different 
values of a^jac- The curves are normalized so that the peak of the ao/ac = 0.01 curve equals 1. 
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